Non-linear effects in the cyclotron resonance of a massless quasi-particle in graphene 
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We consider the classical motion of a massless quasi-particle in a magnetic field and under a weak 
electromagnetic radiation with the frequency w. Due to the non-parabolic, linear energy dispersion, 
the particle responds not only at the frequency u but generates a broad frequency spectrum around 
it. The linewidth of the cyclotron resonance turns out to be very broad even in a perfectly pure 
material which allows one to explain recent experimental data in graphene. It is concluded that the 
linear response theory does not work in graphene in finite magnetic fields. 
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Graphene is a recently discovered two-dimensional (2D) material [US which has attracted much attention due 
to its non-trivial and very interesting physical properties, see reviews [J El- This is a monolayer of carbon atoms 
packed in a dense honeycomb lattice. The band structure of electrons in graphene consists of two bands (electron and 
hole) which touch each other at six corners of the hexagon-shaped Brillouin zone [5] . If graphene is undoped and the 
temperature is zero, the lower (hole) band is fully occupied while the upper (electron) band is empty. The Fermi level 
goes through these six, so called Dirac points. Near these points the states of electron and holes are described by the 
effective Dirac equation with the vanishing effective mass of the quasiparticles, and their energy spectrum is linear, 

£±{p) = ±Vp = ±V\p\. (1) 

Here p = (p x ,Py) is the electronic momentum, counted from the corresponding Dirac points and V 10 8 cm/s is 
the material parameter. It is the massless energy dispersion (JXJ and the Dirac nature of graphene quasiparticles that 
result in its amazing physical properties which attracted so much attention in the past years. 

If a charged massive particle (e.g. electron) is placed in a magnetic field B, it begins to rotate in the plane 
perpendicular to the field B with the frequency uj c = eB/mc, determined by its mass m and charge e. If such a 
rotating particle is irradiated by an external electromagnetic wave with the frequency u>, it absorbs the radiation 
energy if the frequency of the wave is close to the cyclotron frequency, lj ~ oj c . The phenomenon is called the 
cyclotron resonance (CR) and is widely used in solid state physics for characterization of material properties such as 
the charge carrier effective mass and the Fermi surface cross-section in metals and semiconductors [6j . The linewidth 
of the CR absorption line Su) is determined by the scattering rate of 2D electrons and by the radiative decay rate, i.e. 
by losses of energy of the rotating charged particle due to the re-radiation of electromagnetic waves, see e.g. Q. In 
typical semiconductor 2D electron systems, e.g. in GaAs quantum well structures, the linewidth of the CR is small 
as compared to the CR frequency already in magnetic fields of order of 0.1 Tesla (see e.g. Q). 

The magneto-optical conductivity of graphene and its cyclotron response have been studied within the linear 
response theory in a number of publications BBS Emu |13| |. Intensive experimental studies of the CR in graphene 



(as well as other electrodynamic phenomena) have been hampered until recently by the absence of graphene flakes of 
sufficiently large area. With the progress of technology, however, graphene samples of sufficiently large dimensions are 
becoming available now [141 ] , so the growth of experimental activity on the graphene electrodynamics is to be expected. 



First experimental results on the CR in single-layered graphene have been already published in Refs. [15|, llfj (for the 



CR data in bilayer graphene and in thin graphite layers see Refs. [17|, ll8|). In Ref. [15( the infrared spectroscopy has 
been used and in Ref. 11611 t he photoconductive response of graphene has been measured. The characteristic feature 
of both CR experiments |l5l. [l(| is that the CR line was found to be extremely broad. Although the experiments have 
been done in very strong magnetic fields (up to 18 T), the quality factor of the CR line B C /SB C ~lo/5uj only slight! 
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exceeded unity. Following the traditional interpretation of the CR linewidth one had to conclude that in Refs. [151 . 
one dealt with extremely disordered graphene samples. As will be shown below, however, the cyclotron resonance of 
massless quasi-particles ([1]) has very unusual physical properties, and the CR-line can be very broad even in perfectly 
pure graphene. 

The most straightforward way to demonstrate the main idea of this work is to consider the classical motion of a 
quasi-particle with the spectrum {]]) in a uniform magnetic field B = (0, 0, B) and in the presence of an external 
electromagnetic wave. The equation of motion reads 

I— ivxB-dKt), v = v£. (2) 
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FIG. 1: (Color online) Time dependencies of the momentum of the massive (a) and the massless (b) quasi-particles under the 
CR conditions. 



If the external electric field E(t) is zero, the particle rotates in the magnetic field with the energy-dependent cyclotron 
frequency lj c — eBV/pc = eBV 2 /c£, where the energy £ = Vp is the integral of motion (e.g. [19| . §21). If the 
external electric field is switched on, one usually assumes that the particle absorbs the electromagnetic wave energy 
at the cyclotron frequency u> c . This is perfectly true for electrons with the parabolic spectrum, but for the massless 
graphene quasi-particles |T]) the situation turns out to be much more complicated and more interesting. 

In Figure Q] we show results of the numerical solution of Eq. © with the initial conditions p|t=o = Po- We 
assume that the radiation frequency uj coincides with the cyclotron frequency lo c and that the external electric field 
Eo(i) = Eo(t)(cosojt, sinut) is circularly polarized in the direction coinciding with the sense of the cyclotron rotation of 
the particle. The spectrum of the incident radiation Si nc (Q) thus consists of only one spectral line, Si nc (fl) cx 5(Q— uS), 
and the particle is under the perfect CR conditions. For the sake of clarity we further neglect all dissipative processes 
such as the scattering of electrons by impurities and phonons, as well as the radiative decay ([Uj], §67). We also 
assume that the external electric field Eq is weak, T = eEo/tupo <C 1, meaning that the energy absorbed by the 
particle during one oscillation period is small as compared to its average energy. Usually this corresponds to the 
linear response regime. 

Figure [T^, shows the known results for a standard (massive) particle (in this case v = p/m in (|2|) and lu c — eB/mc). 
Under the chosen conditions, the massive particle is always in resonance with the external field, it continuously absorbs 
the radiation energy, and the absolute value of its momentum p linearly grows with time. Rotating with the frequency 
u> the massive particle re-emits (scatters) the radiation of the incident wave at the same frequency u>, and the spectrum 
of the scattered waves S sca t(Q) is also proportional to 6(Q — oj). 

The motion of the massless particle is essentially different. Its cyclotron frequency uo c — eBV/pc = eBV 2 /£c 
depends on its energy. Initially, when the external field is switched on, the particle is also in resonance with the 
external radiation, lo = lo c — eBV/poc, and starts to get energy from the wave. But, as soon as its energy increases, 
p > po, it gets out of the resonance and ceases to absorb the radiation energy. Its energy £ = Vp and the absolute 
value of the momentum p decreases and the particle again turns out to be in the resonance with the wave. Then the 
process repeats itself, Fig. [TJj, and the particle energy oscillates in time with the period depending on the external 
field amplitude Eq- The time dependence of the momentum of the massless particle is not a harmonic function with 
only one Fourier component fi = u>. It also contains other (higher and lower) harmonics. 

To study this phenomenon further we introduce, instead of p x (t) and p y (t), the new variables p(t) and <f>(t) according 
to the formulas 

( - n u\ ( -M"t + 4>(t)) 

\Py(t) J ~ P[ > \ cosM + ^i)) 
and investigate the time dependencies of the momentum p(t) and the phase 4>{t). Equations ^ are then rewritten as 

p(t) = eE sin (t)(t), (4) 




p(t)<j)(t) = -ujp(t) + eVB/c + eE Q cos </>(£), 



(5) 
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and the initial conditions are p\t=a — Pa, <j>\t=o = 4>a- Figures [2^,b show the time dependencies of p(t) and 4>(t) at 
different values of the electric field parameter T = eEo/ojpo- Both functions are modulated, and the modulation 
frequency, as well as the amplitude of the momentum oscillations, decrease with T . In contrast, the amplitude of the 
phase oscillations remains independent on the electric field Eq, Fig. [2b- 

The behavior of the p(t) and 4>{t) oscillations also depends on the initial phase (f>o, Fig. [2j;,d. One sees that not 
only the oscillation period and amplitude may depend on <f>Q, but even the overall shape of oscillations, see the case 
4>0 = 7r (green dash-dotted curve on Fig. [2j:,d). The linear contribution to the phase, Fig. [2li, means the frequency 
shift, since <j>(t) = —at + <f)(t) leads to Lot + cf){t) = (w — a)t + <f>(t). 

So far we discussed the behavior of the momentum of the massless quasi-particle. The experimentally measured 
value, however, is the velocity 

W p(t) { co S (u>t + ^t)) )■ w 

Its absolute value |v(t)| = V remains constant independent on how big or small the particle momentum is. This 
means that, even in a weak external electric field Eq which produces small momentum oscillations, the induced ac 
current j(t) = — en s v(t) (in a system of many massless quasi-particles) can be rather strong. For example for a 
typical electron density of n s — 10 12 cm~ 2 in graphene the estimate gives a very big value of the current j m 16 
A/cm. Another important feature of the time dependent velocity ((6J is that its phase remains to be modulated even 
in the limit of week electric fields F — > 0, Fig. Od. The Fourier spectrum of the velocity, Figure shows therefore a 
broad spectrum around the central frequency Q = lo at T = eEo/wpo <C 1. Since the velocity of particles determines 
the current and the amplitudes of the scattered electromagnetic waves, Fig. [3] shows, in fact, the Fourier spectrum 
of the scattered waves S sca t(Q). The function S sca t(£l) contains a lot of Fourier harmonics although the spectrum 
of the incident wave is proportional to 5(Q — lo). We emphasize that we discuss the CR resonance of the massless 
quasi-particle in the absence of any scattering. 

In the limit T -C 1 one can get some more results analytically. Substituting pit) =po[l + with \q(t) \ < 1, we 
rewrite Eqs. Q-© as 

q(t) = LoFsincj), (f)(t) = — uiq + uj J 7 cos (j). (7) 

It will be seen from the result © that q oc \jlP , therefore at T 1 one can neglect the last term in the second 
equation ([7]) as compared to the first one. Then the problem is reduced to the equation of the non-linear pendulum 

0(f) = -u?T sin <j>. (8) 

If the initial phase </>o is small, 4>q ^ 1, the solution assumes the form 



<f>(t) — (j>o cos(u>Vrt) , q(t) — (f>o\^sm(uiVlFt). 



(9) 
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FIG. 3: Frequency spectra of the velocity v x for T = eEo/ujpo = 0.01 and 4>o = t/2 (a) and 4>o = n (b). 



The amplitude and the phase of the momentum oscillations are thus modulated with the frequency uj\[T proportional 
to the square root of the electric field. 

Substituting the phase <j>(t) from d9j) to (J6j) one can calculate the Fourier spectrum of the velocity v(fi) = 
J \-(t)e~ int dt/2n. For example, for the v y component we get 

oo oo 

v y (n) =2-2] ik MM$(n - w - kuVr) + — {-i) k Jk(MS(n + w- kuVf). (10) 

fe=— oc h— — oo 

The spectrum consists of an in infinite number of satellite harmonics at the frequencies f2 = uifcwVf with k = 1, 2, . . .. 
The amplitudes of these harmonics Jk(4>o) are determined only by the initial phase of the particle and do not depend 
on the electric field. 

So far we discussed the response of only one quasi-particle and found that it is complicated and very different even 
for particles with the same po but different initial phases 4>o- In a real graphene system one deals with many particles 
with different initial phases and initial momenta po <pf lying inside the Fermi circle. The electromagnetic response 
of the whole many-particle system will be a superposition of partial responses of all the individual particles and will 
therefore be very complicated. The time dependence of the current (see example in Figure 2J looks very chaotically, 
and the scattered wave spectrum S sca t(Q) will have a broad peak around Q = lu with the quality factor only slightly 
exceeding unity. Such a large linewidth of the scattered wave resonance is not related to the scattering processes but is 
due to the linear graphene dispersion (Q]) . Notice that the very broad CR lines have been observed in the experiments 
TEl 16j |. As follows from our results they are not necessarily related to disorder: The broad CR line can be observed 
even in a perfectly pure graphene. It should also be noticed that chaotic motion of graphene quasi-particles has been 
discussed in other respects in Refs. (20L I21I and the non-linear features of the CR in systems with the non-parabolic 
energy dispersion have been studied in [22j. 

Above we assumed that the external ac electric field is weak, T <C 1. In such cases one usually uses the linear 
response theory which automatically implies that the system responds at the frequency coinciding with the frequency 
of the external radiation. As we have shown, this is not the case for the massless quasi-particles with the spectrum 
$Q . It should therefore be concluded that the current-vs-field dependence in graphene is non-analytical at low electric 
fields, and that the linear response theory does not adequately describe the electromagnetic response of graphene in 
finite magnetic fields. 

Mathematically, the non-analyticity of the graphene response function originates from the singular Lorentz-force 
term in the equation of motion ([2]). This singularity disappears at B — 0, therefore in zero magnetic field the linear- 
response theory should be valid. This agrees with our previous results from Refs. [23L I24] ] where the non-linear 
frequency multiplication effect, predicted at T S> 1, was shown to disappear in the weak fields. 

To summarize, electromagnetic response of graphene in a magnetic field was shown to be essentially non-linear. The 
system irradiated by the monochromatic wave scatters the radiation in a broad frequency range. In agreement with 
experiments, the CR line was shown to be very broad in graphene, even if the intensity of the incident radiation is 
weak and all scattering processes are neglected. The linear response theory was shown to be not completely adequate 
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time, atV2n 

FIG. 4: (Color online) Time dependence of the current in a system of 48 massless particles. It is assumed that the particles are 
Fermi-distributed in the p-space at t = and that the electric field is zero at iot/2n < 5 and smoothly grows up to a constant 
value at 5 < ujt/2-n < 10. 

for the description of electromagnetic properties of graphene. Instead, the methods of the non-linear dynamics and 
the theory of chaos should be used for this purpose. 
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